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A 2 

We investigate the occurrence of power terms ^ in the running QCD coupling by analysing non-perturbative 
measurements of Q. a (p) at quite low momenta obtained from the lattice three-gluon vertex. Our study provides 
some evidence for such a contribution. The phenomenological implications of such a presence are reviewed. 
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1. Introduction and motivations 

The Operator Product Expansion (OPE) is the 
standard way to parametrise non-perturbative 
QCD effects in terms of power corrections, the 
powers involved beeing uniquely fixed by the op- 
erator content of the theory. Anyway, due to the 
asymptotic nature of QCD perturbative expan- 
sions, i.e. due to renormalons ambiguities, power 
corrections are reshuffled between operators and 
coefficient functions |3J . It has recently been con- 
jectured Q that power corrections which are not 
a priori expected from OPE may appear in the 
expansion of physical observables, via much the 
same analysis as for renormalons, if one allows the 
presence of (UV-subleading) power corrections to 
ct s (j>). The lattice community knows that per- 
turbative logarithms are of course not the only 
contribution to the running coupling. Aim of 
the present work is to search for non-perturbative 
contributions to the running coupling in the form 
of power corrections, with the given power (^r). 
Before briefly reviewing some theoretical argu- 
ments in support of this theoretical prejudice, 
a couple of considerations are in order. While 
power corrections to a s (p) arise naturally in many 
physical schemes HQ and their occurrence can- 
not be excluded a priori in any scheme, their 
non-perturbative nature makes it very hard to as- 



sess their scheme dependence, which is only very 
weakly constrained by the general properties of 
the theory. 

1.1. A lesson from the Gluon Condensate 

A stage on which all the considerations above 
apply is provided by the efforts towards a lattice 
determination of the Gluon Condensate, which is 
given in terms of Wilson loops W J|] . The OPE 
for W is to be written down as 



W = Wo + 



A 4 W 4 

Q 4 



(1) 



Wo beeing the contribution related to the identity 
operator and the second term beeing the "gen- 
uine" (dim=4) condensate (Q is the scale; on the 
lattice Q ~ 1/a). Perturbative contributions are 
present only in Wo, so that a standard procedure 
to extract the condensate was to subtract from 
MonteCarlo measurements of Wilson loops their 
perturbative expansions. From general consider- 
ations the expected form of Wq is 



Wo 



' dp 2 p 2 . 

pA 2 P H 



2 a s ( 
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A renormalon analysis of this formula readily 
shows that a 0(A 4 /Q 4 ) ambiguity is present in 
Wq , thus preventing any unambiguous result from 
the above procedure. The situation is anyway 
even more intricate: actually, after having re- 
summed Wq (within the mentioned ambiguity), 
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performing the subtraction leaves one with some- 
thing that scales like 0(A 2 /Q 2 ) |§- The puzzle 
can be sorted out if one supposes in (||) contribu- 
tions of the form a s (p) ~ • 

1.2. Static quarks potential and confine- 
ment 

a 2 

A nice argument to support the presence of 
contributions to a s (p) comes from considerations 
involving confinement. Consider the interaction 
of two heavy quarks in the static limit || . Within 
the Born approximation one can obtain the static 
potential from 



V(r) 



oc 
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The above formula has been written down in 
the renormalon-style: plugging in a constant a s 
yields the Coulomb potential (V(r) w 1/?"), while 
the leading-logs expression for a s generates vari- 
ous power corrections to the potential; however 
the string tension contribution V(r) ~ Kr is 
missing. Such a contribution is obtained by plug- 
ging in a power correction to a s of the form we 
are looking for. Notice that the above consid- 
erations are in a sense not new in the context 
of non-perturbative contributions to the running 
coupling Consider the "force" definition of 
the running coupling: 



(4) 



where again V(r) is the static interquark poten- 
tial. By keeping into account the string tension 
contribution to V(r), one obtains a 1/Q 2 contri- 
bution, whose order of magnitude is given by the 
string tension itself. While this term was mainly 
considered as a sort of ambiguity, resulting in an 
indetermination in the value of a(Q) at a given 
scale, it can be interpreted as a clue for the ex- 
istence of a contribution, providing also an 
estimate for the expected order of magnitude of 
it, at least in one (physically sound) scheme. 

2. Fits to lattice data 

The lattice data for a s (p) that we used for 
our analysis were obtained by evaluating two- 



and three-point off-shell Green's functions of the 
gluon field in the Landau gauge and imposing 
non-perturbative renormalisation conditions on 
them, for different values of the external mo- 
menta ||. Such a definition of the coupling 
corresponds to a momentum-subtraction renor- 
malisation scheme in continuum QCD |9). The 
data available at the moment are quite noisy 
and lattice artifacts are still to be fully assessed 
(they seem under control from the analysis in j8|). 
Given the particular ansatz to which we want to 
fit the data (perturbative expressions of given or- 
der plus a power correction), a peculiar momen- 
tum interval has to be singled out. On one hand, 
the momentum range should start well above the 
location of the perturbative Landau pole, but it 
should nonetheless include low scales where power 
corrections may still be sizeable. On the other 
hand, the requirement of keeping the effects of 
the finite lattice spacing under control in the nu- 
merical data for a s induces a natural UV cutoff. 

2.1. Two-loop analysis 

At two-loop level we fit data to the formula 



I 
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h log(L 2 ) 
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A 2 
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where L 2 = log^/A 2 .;). Notice that since the 
value of A is expected to carry a sizeable depen- 
dence on the order of the perturbative calcula- 
tion, we append a subscript. We were able to 
check that out of power corrections of the form 



{~^~Y the value z 



1 is indeed the best choice, 
which thing we interpret as a confirmation of our 
theoretical prejudice. Our best fit singles out 
the values (A 2/ = 0.84(l),c 2 ; = 0.31(3)), with a 
Xdof 5: 1-8. The momentum interval which is best 
described is p ~ 2 — 3 GeV, fully consistent with 
the requirements we have already mentioned. 

2.2. Three-loop analysis 

Given the possible interplay between a descrip- 
tion in terms of power corrections and our ig- 
norance about higher orders, a three-loop anal- 
ysis is compelling. A major obstacle for such an 
analysis is actually the fact that the first non- 
universal /3-function coefficient b 2 is not known 
for our scheme. Therefore we also fitted data to 
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the formula 



a s (p) 



bo L 3 



h \og(L 3 ) 1 



' h eff h 2 ^ 

^- + ^(log 2 (L 3 )~log(L 3 ) + l) 
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A 2 



(6) 



where L3 = log(p 2 /A'^ l ) and b 2 is an effec- 
tive /3-function coefficient to be determined from 
the fit. In order to gain insight, we started by 
putting C31 — 0, obtaining the values (A3; = 
0.72{l),b e 2 ff = 1.3(1)), with x 2 do} ~ 1-8, in a mo- 
mentum range p ~ 2 — 3 GeV(dashed line in the 
figure). The value of is expected to be a 
reliable estimate of 62 , as confirmed by simple ar- 
guments concerning the convergence properties of 
the expansion of our coupling in powers of other 
couplings. What is more important, the value of 
A3/ agrees with other determinations via pertuba- 
tive matching j[o|. By keeping into account also 
C3/, moving around the same range for A3; one 
obtains the values (A 3; = 0.72(1), b e 2 ff = 1.0(1), 
C3Z = 0.41(2)), yielding a Xdof ~ 1-8 in a momen- 
tum range p ~ 1.8 — 3 GeV(solid line in the fig- 
ure). While &2 s tiU makes sense, it emerges that 
c 2 /A 2 , = 0.22(2) GeV 2 ~ c 3 zA 2 , = 0.21(2) GeV 2 , 
which is just the order of the standard estimate 
for the string tension. 



3. Conclusions and perspectives 

Some preliminary evidence in support of a ^ 
contribution to a s (p) has been found, which ap- 
pears at this stage disentangled from possible 
perturbative ambiguities: perturbative and non- 
perturbative (power) contributions do not mix in 
our formulae when upgrading from a two-loop to 
a three-loop description. While our results need 
to be further tested, still they provide evidence in 
support of the conjecture that "anomalous" (i.e., 
not accounted for by OPE) power corrections into 
current correlation functions and physical observ- 
ables may be generated by power corrections to 
a s (p). 




Figure 1. Three-loop fits to lattice data for the 
coupling, with and without power corrections. 
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